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Abstract 



Dixon's multipoles for a system of N relativistic positive-energy scalar par- 
ticles are evaluated in the rest-frame instant form of dynamics. The Wigner 
hyperplanes (intrinsic rest frame of the isolated system) turn out to be the 
natural framework for describing multipole kinematics. In particular, con- 
cepts like the barycentric tensor of inertia can be defined in special relativity 
only by means of the quadrupole moments of the isolated system. 
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I. INTRODUCTION. 



In a recent paper []I| we have given a complete treatment of the kinematics of the rela- 
tivistic N-body problem in the rest-frame instant form of dynamics We have shown 

in particular how to perform the separation of the center-of-mass motion in the relativistic 
case. This requires the reformulation of the theory of relativistic isolated systems on arbi- 
trary spacelike hypersurfaces Q. This description is also able to incorporate the coupling to 
the gravitational field. This is essentially Dirac's reformulation P of classical field theory 
(suitably extended to particles) on arbitrary spacelike hypersurfaces {equal time surfaces) 
and provides the classical basis of the Tomonaga-Schwinger formulation of quantum field 
theory. For each isolated system (containing any combination of particles, strings and fields) 
one gets a reformulation as a parametrized Minkowski theory [0, with the extra bonus of 
having the theory already predisposed to the coupling to gravity in its ADM formulation, 
but with the price that the functions z^(r, a) describing the embedding of the spacehke 
hypersurface in Minkowski spacetime become configuration variables in the action principle. 
Since the action is invariant under separate r-reparametrizations and space-diffeomorphisms, 
there emerge first class constraints ensuring the independence of the description from the 
choice of the 3+1 splitting. The embedding configuration variables z^{T^a) are the gauge 
variables associated with this kind of general covariance. 

Let us remark that, since the intersection of a timelike worldline with a spacelike hyper- 
surface corresponding to a value r of the time parameter is identified by 3 numbers a = f]{T) 
instead of four, in parametrized Minkowski theories each particle must have a well defined 
sign of the energy: therefore we cannot describe the two topologically disjoint branches of 
the mass hyperboloid simultaneously like in the standard manifestly Lorentz-covariant the- 
ory. As a consequence there are no more mass-shell constraints. Therefore, each particle 
with a definite sign of the energy is described by the canonical coordinates rfj(r), Kj(r) with 
the derived 4-position of the particles given by (r) = 2;'^(r, fjiij)). The derived 4-momenta 
(r) are /5j-dependent solutions of pj — emf = with the chosen sign of the energy. 

As said, parametrized Minkowski theories have separate spatial and time reparametriza- 
tion invariances, which imply the independence of the description from the choice of the 
3+1 splitting of Minkowski spacetime. In Minkowski spacetime we can restrict the foliation 
to have spacelike hyperplanes as leaves. Then, for each configuration of the isolated system 
with timelike 4-momentum, we can restrict ourselves to the special foliation whose leaves are 
the hyperplanes orthogonal to the conserved system 4-momentum, which have been named 
Wigner hyperplanes. This special foliation is intrinsically determined only by the configu- 
ration of the isolated system. In this way it is possible to arrive at the definition of the 
Wigner- covariant rest-frame instant form of dynamics for every isolated system whose con- 
figurations have well defined and finite Poincare generators with timelike total 4-momentum 
(see Ref. for the traditional forms of dynamics). 

This formulation clarifies the non-trivial definition of a relativistic center of mass. As 
well known, no such definition can enjoy all the properties of the non-relativistic center of 



^Such hypersurfaces are the leaves of a foliation of Minkowski spacetime (namely one among its 
3+1 splittings) and are equivalent to a congruence of timelike accelerated observers. 
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mass. See Refs. [|,|],^T]-|T^ for a partial bibliography of all the existing attempts. 

As shown in Appendix A of Ref. |l| only four first class constraints survive in the rest- 
frame instant form on Wigner hyperplanes and the original configuration variables -2''(t, a), 
fjiir) and their conjugate momenta p^{T,a), Riij) are reduced to: 

i) a decoupled particle 5;^(r), (the only remnant of the spacelike hypersurface) with a 
positive mass = y^epf determined by the first class constraint eg — Mgys ~ and with its 
rest-frame Lorentz scalar time T, = put equal to the mathematical time as the gauge 
fixing Ts — r ^ to the previous constraint. Here, x'^{t) is a non-covariant canonical vanahle 
for the 4-center of mass. After the elimination of T, and es with the previous pair of second 
class constraints, one remains with a decoupled free point {point particle clock) of mass Mgys 
and canonical 3-coordinates Zg = €s{xs — ^x°), ks = The non-covariant canonical a;^(r) 
must not be confused with the 4- vector x^(r) = z^{T,(f = 0) identifying the origin of the 
3-coordinates a inside the Wigner hyperplanes. The worldline a;^(r) is arbitrary because 
it depends on x^(0) and its 4- velocity x^(r) depends on the Dirac multipliers associated 
with the 4 left first class constraints 0, as it will be shown in the next Section. The unit 
timelike 4- vector u^{ps) = P^/^s is orthogonal to the Wigner hyperplanes and describes their 
orientation in the chosen inertial frame. 

ii) the particle canonical variables ffiir), Ri{T) [| inside the Wigner hyperplanes. They are 
restricted by the three first class constraints (the rest-frame conditions) /?+ = Y^f=i Ki ~ 0. 
Since the role of the relativistic decoupled ^-center of mass is taken by x^{t) (or, after the 
gauge fixing — r 0, by an external 3-center of mass Zs, defined in terms of and 
|l|]), the rest-frame conditions imply that the internal 3-center of mass g+ = acom is a gauge 
variable, which can be eliminated with gauge fixings 0. 

Therefore, we need a doubling of the concepts: 

1) There is the externa/ viewpoint of an arbitrary inertial Lorentz observer, who describes 
the Wigner hyperplanes, leaves of a foliation of Minkowski spacetime, determined by the 
timelike configurations of the isolated system. A change of inertial observer by means of a 
Lorentz transformation rotates the Wigner hyperplanes and induces a Wigner rotation of 
the 3-vectors inside each Wigner hyperplane. Every such hyperplane inherits an induced 
internal Euclidean structure while an external realization of the Poincare group induces the 
internal Euclidean action. As said above, an arbitrary worldline [centroid) a;^(r) is chosen 



'^Msys being the invariant mass of the isolated system. 

'^Zg/es is the classical analogue of the Newton- Wigner 3-position operator, only covariant under 
the Euclidean subgroup of the Poincare group. 

^Therefore this arbitrary worldline may be considered as an arbitrary centroid for the isolated 
system. 

^They are Wigner spin-1 3-vectors, like the coordinates a. 

^For instance g+ ~ implies that the internal 3-center of mass is put in the origin x^(r) = 
z^'{T,a = 0). 
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as origin of the internal 3-coordinates on the Wigner hyperplanes. 

Three external concepts of 4-center of mass can be defined by using the external reahza- 
tion of the Poincare algebra (to each one there corresponds a 3-location inside the Wigner 
hyperplanes) : 

a) the external non-covariant canonical ^-center of mass (also named ^--center of spin) 
(with 3-location a), 

b) the externa/ non-covariant non-canonical M0ller ^-center of energy (with 3-location 

c) the external covariant non-canonical Fokker-Pryce 4-center of inertia F/* (with 3- 
location cry). 

Only the canonical non-covariant center of mass x^(r) is relevant in the Hamiltonian 
treatment with Dirac constraints, while only the Fokker-Pryce is a 4-vector by con- 
struction. See Ref. [||] for the construction of the ^--centers starting from the corresponding 
3-centers (3-center of spin |T^], 3-center of energy |]TU||, 3-center of inertia [IT],] 



2) There is the internal viewpoint inside the Wigner hyperplanes associated to a un- 
faithful internal realization of the Poincare algebra: the total internal 3-momentum of the 
isolated system vanishes due to the rest-frame conditions. The internal energy and angular 
momentum are the invariant mass Mgyg and the spin (the angular momentum with respect 
to x^(t)) of the isolated system respectively. With the internal realization of the Poincare 
algebra we can define three internal 3-centers of mass: the internal canonical 3-center of 
mass, the internal M0ller 3-center of energy and the internal Fokker-Pryce 3-center of in- 
ertia. But, due to the rest-frame conditions, they coincide. As a natural gauge fixing to 
the rest-frame conditions we can add the vanishing of the internal Lorentz boosts: it is 
equivalent to locate the internal canonical 3-center of mass (f+ in c? = 0, i.e. in the origin 
x^(r) = z^{t, 0). With these gauge fixings and with Tg — r^O, the worldline x^(r) becomes 
uniquely determined except for the arbitrariness in the choice of x^(0) [u^{ps) = Ps/^s] 

x^,{t)=x^{0) + u^{Ps)T,, (1.1) 

and coincides with the external covariant non-canonical Fokker-Pryce 4-center of inertia, 
x^{r)=x^iO) + Y^^ i- 

This doubling of concepts with the external non-covariant canonical 4-center of mass 
x^(r) (or with the external 3-center of mass Zg when Ts — t ^ 0) and with the internal 
canonical 3-center of mass g+ ~ replaces the separation of the non-relativistic 3-center of 
mass due to the Abelian translation symmetry. The non-relativistic conserved 3-momentum 
is replaced by the external Ps = egkg, while the internal 3-momentum vanishes, /?+ ~ 0, as 
a definition of rest frame. 

In the final gauge we have = Mgys, Tg = t and the canonical basis i^, /cg, i/j, Ki is 
restricted by the three pairs of second class constraints k+ = Y^iLi ~ 0, g+ ~ 0, so that 
6N canonical variables describe the N particles like in the non-relativistic case. We still need 
a canonical transformation fji, Ki \—>- 0], 0], Pa, vr^ [a = 1, .., — 1] identifying 

a set of relative canonical variables. The final 6N-dimensional canonical basis is Zg, kg, pa, 
TTa- To get this result we need a highly non-linear canonical transformation |]l[, which can 
be obtained by exploiting the Gartenhaus-Schwartz singular transformation pi|. 



In the end we obtain the Hamiltonian for relative motions as a sum of N square roots, 
each one containing a squared mass and a quadratic form in the relative momenta, which 
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goes into the non-relativistic Hamiltonian for relative motions in the hmit c ^ oo. This 
fact has the following implications: 

a) if one tries to make the inverse Legendre transformation to find the associated La- 
grangian, it turns out that, due to the presence of square roots, the Lagrangian is a hy- 
perelliptic function of already in the free case. A closed form exists only for N=2, 

mi = m2 = m: L = —em\J 4 — . This exceptional case already shows that the existence of 
the limiting velocity c (or in other terms the Lorentz signature of spacetime) forbids a linear 
relation between the spin (center-of-mass angular momentum) and the angular velocity. 

b) the N quadratic forms in the relative momenta appearing in the relative Hamiltonian 
cannot be simultaneously diagonalized. In any case the Hamiltonian is a sum of square 
roots, so that concepts like reduced masses, Jacobi normal relative coordinates and tensor of 
inertia cannot be extended to special relativity. As a consequence, for example, a relativistic 
static orientation-shape SO (3) principal bundle approach can be implemented only by 
using non- Jacobi relative coordinates. 

c) the best way of studying rotational kinematics (the non-Abelian rotational symmetry, 
associated with the conserved internal spin) is based on the canonical spin bases with the 
associated concepts of spin frames and dynamical body frames introduced in Ref. they 



can be build in the same way as in the non-relativistic case |]I6[ starting from the canonical 

basis Pa, Tla- 

Let us clarify this point. 

In the non-relativistic N-body problem it is easy to make the separation of the absolute 
translational motion of the center of mass from the relative motions, due to the Abelian na- 
ture of the translation symmetry group. This implies that the associated Noether constants 
of motion (the conserved total 3-momentum) are in involution, so that the center-of-mass 
degrees of freedom decouple. Moreover, the fact that the non-relativistic kinetic energy of 
the relative motions is a quadratic form in the relative velocities allows the introduction of 
special sets of relative coordinates, the Jacobi normal relative coordinates, which diagonalize 
the quadratic form and correspond to different patterns of clustering of the centers of mass 
of the particles. Each set of Jacobi normal relative coordinates organizes the N particles 
into a hierarchy of clusters, in which each cluster of two or more particles has a mass given 
by an eigenvalue {reduced masses) of the quadratic form; Jacobi normal coordinates join the 
centers of mass of pairs of clusters. 

However, the non-Abelian nature of the rotation symmetry group, whose associated 
Noether constants of motion (the conserved total angular momentum) are not in involution, 
prevents the possibility of a global separation of absolute rotations from the relative motions, 
so that there is no global definition of absolute vibrations. This has the consequence that an 
isolated deformable body can undergo rotations by changing its own shape (see the examples 
of the falling cat and of the diver). It was just to deal with these problems that the theory of 
the orient at ion- shape SO (3) principal bundle approach has been developed. Its essential 
content is that any static (i.e. velocity-independent) definition of body frame for a deformable 



^See Ref. Q for a review of this approach used in molecular physics for the definition and study 
of the vibrations of molecules. 
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body must be interpreted as a gauge fixing in the context of a SO (3) gauge tlieory. Botfi 
tlie laboratory and tlie body frame angular velocities as well as the orientational variables 
of the static body frame become thereby unobservable gauge variables. This approach is 
associated with a set of point canonical transformations, which allow to define the body 
frame components of relative motions in a velocity-independent way. 

Since in many physical applications (e.g. nuclear physics, rotating stars,...) angular 
velocities are viewed as measurable quantities, it is desiderable to have an alternative for- 
mulation complying with this requirement, possibly generalizable to special relativity. This 
program has been realized in a previous paper |TB], of which Ref. is the relativistic 



extension. Let us summarize the main points of our formulation. 

First of all for > 3 we have constructed in Ref. [|l6l a class of non-point canonical 



transformations, which allows to build the canonical spin bases quoted above, which are 
connected to the patterns of the possible clusterings of the spins associated with relative 
motions. The definition of these spin bases is independent of Jacobi normal relative coordi- 
nates, just as the patterns of spin clustering are independent of the patterns of center-of-mass 
Jacobi clustering. We have found two basic frames associated to each spin basis: the spin 
frame and the dynamical body frame. Their construction is guaranteed by the fact that, 
besides the existence on the relative phase space of a Hamiltonian symmetry left action of 
SO(3)00, it is possible to define as many Hamiltonian non-symmetry right actions of S0(3) 
as the possible patterns of spin clustering. While for N=3 the unique canonical spin basis 
coincides with a special class of global cross sections of the trivial orientation-shape SO (3) 
principal bundle, for > 4 the existing spin bases and dynamical body frames turn out 
to be unrelated to the local cross sections of the static non-trivial orient at ion- shape SO (3) 
principal bundle, and evolve in a dynamical way dictated by the equations of motion. In this 
new formulation both the orientation variables and the angular velocities become measurable 
quantities in each canonical spin basis by construction. 

For each N every allowed spin basis provides a physically well-defined separation between 
rotational and vibrational degrees of freedom. The non-Abelian nature of the rotational sym- 
metry implies that there is no unique separation of absolute rotations and relative motions. 
The unique body frame of rigid bodies is replaced by a discrete number of evolving dynami- 
cal body frames and of spin canonical bases, both of which are grounded on patterns of spin 
couplings, direct analogues of the coupling of quantum angular momenta. 

In this paper we complete our study of relativistic kinematics for the N-body system by 
evaluating its rest-frame Dixon multipoles ||17||p^. Let us recall that this method is being 
used for treating extended systems in astrophysics. The starting point is the definition of 



^We adhere to the definitions used in Ref. |15]; in the mathematical literature our left action is a 
right action. 

^Their generators are the center-of-mass angular momentum Noether constants of motion. 
^^Their generators are not constants of motion. 

^^See Ref. [18| for the definition of Dixon's multipoles in general relativity. 
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the energy momentum tensor of the N positive-energy particles on the Wigner hyperplane. 
In this way we see that the Wigner hyperplane is the natural framework for reorganizing 
a lot of kinematics connected with multipoles. Moreover, only in this way, a concept like 
the barycentric tensor of inertia can be introduced in special relativity by means of the 
quadrupole moments. 

A review of the rest-frame instant form of dynamics for N scalar free positive-energy 
particles is given in Section II. 

In Section III we evaluate the energy momentum tensor on the Wigner hyperplanes. 

Dixon's multipoles for the system are defined in Section IV. A special study of monopole, 
dipole and quadrupole moments is given together with the multipolar expansion. 

Other properties of Dixon's multipoles are reviewed in Section V. 

Some comments on open problems are given in the Conclusions. 

The non-relativistic N-particle multipolar expansion is given in Appendix A, while in 
Appendix B there is a review of symmetric trace-free (STF) tensors. 
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II. REVIEW OF THE REST-FRAME INSTANT FORM. 



Let us review the system of N free scalar positive-energy particles in the framework of 
parametrized Minkowski theory (see Appendices A and B and Section II of Ref. As 
said in the Introduction, each particle is described by a configuration 3- vector fjiij). The 
particle worldline is x^{t) = z'^{T,r]i{T)), where z^{T,a) are the embedding configuration 
variables descibing the spacelike hypersurface S,- 0. 

The system is described by the action 



S = J dTdPaC{T,a) = J dTL{T), 



N 



i=l 

N 



L{r) = -Y,mi^g^Jj~^~{j^ (2.1) 



i=l 

The action is invariant under separate r- and a-reparametrizations. 
The canonical momenta are 



Kif ( T I 



dLjr) 

^rr(r, fjiir)) + gfsjT, r]i{T))r]f{T) 

\/gTAr, f]i{T)) + 2grf{T, f]t{T))r]l{T) + gfsir, Vi{r))r]l {r)r]! [t) 



rrii 



{n!{T),^,,(T)} = ~6,,St. (2.2) 
The canonical Hamiltonian He is zero, but there are the primary first class constraints 



^^The foliation is defined by an embedding R x T, ^ M^, ij-,'^) ^ -2^' (t, i?) , with S an abstract 
3-surface diffeomorphic to B?. S,- is the Cauchy surface of equal time. The metric induced on it is 
QabIA — ^A^nuz'^ on S,-, a functional of z^, and the embedding coordinates z^{t, a) are considered 
as independent fields. We use the notation = (t,^^ of Ref s. |J|. The z^(a) = dz''{a)/da^ 
are flat cotetrad fields on Minkowski spacetime with the z^'s tangent to S,-. While in Ref. Q we 

used the metric convention -q^i, = e(H ) with e = ±, in this paper we shall use e = 1 like in 

Ref. §. 
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^ I 

n^{T,a) = p^,{T,a) - l^{T,a)Y,^^i.^ - mi.'r))^ T^f - (r, ct) (r ) /t^^ (r ) - 

i=l 
N 

- Zr^ir, a)Y'ir, a) - U'r)Ws ~ 0, (2.3) 

i=l 

SO that the Dirac Hamiltonian is Ho = J d^aX'^{T, a)T-Cf^{T, a), with the A'^(r, a) Dirac mul- 
tiphers. 

The conserved Poincare generators are (the suffix "s" denotes the hypersurface E,-) 



Jr = / d'a[z^{T, a)p''{T, a) - z'^{t, a)p^{T, a)]. (2.4) 

After the restriction to spacehke hyperplanes the Dirac Hamiltonian is reduced to = 
A^(r)7Y'^(r) + A^i,(r)7Y^'^(r) (only ten Dirac multipliers survive) with the remaining ten 
constraints given by 

N N 



i=l i=l 

W''{t) = h'^.{T) J d^a(T'n''{T,a) - b;{T) J d^^aa^n^ir, a) = 



N 



i=l 

N 

[6^(r)6Kr)-6,^(r)6^(r)]E^rM<(r)^0. (2.5) 



i=l 



Here 5^*^ is the spin part of the Lorentz generators 

5f = 6^(r) j d^'aa'p-'iT, a) - 6^(r) J d^aa'p^'ir, a). (2.6) 
On the Wigner hyperplane^^ we have the following constraints and Dirac Hamiltonian 

mn 

N N 

n^{T) =p^- u^{p,) E Vmf + Rf + e^{u{ps)) E ^ir = 

i=l i=l 



^^On it we use the notation A = {T,r). The 3- vectors B = {b^} on the Wigner hyperplanes are 
Wigner spin-1 3- vectors. 

^^e(f(n(ps)) = L^^r{Ps,Ps) s-^d ^ri'^iPs)) = u^{Ps) = L^o{Ps,Ps) ^-^e the columns of the standard 
Wigner boost for timelike Poincare orbits. See Appendix B of Ref. 111. 
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N N 

4=1 i=l 

or 

^ I 

- Msys ~ 0, Msys = TJlf + 

N 

Psys = ^+ = ~ 0, 

i=l 

Ho = X^{T)n,{T) = X{T)[e, - 

A(r) ^ -x,^(r)M^(p,), 
A,.(r) ^ -x,^,(r)e^(M(p,)), 

C(r) = -A(r)«'^(p,), 

x^{t) = {x^{t),Ho} = Ur){x^{r),mr)} ^ 

-A^(r) = -X{T)u^{p.,) + e^{u{ps))Xr{r), 

X^{t)~X\t)>0, = -A(r), 

i^jr) ^ -A(r)n^(p,) + A,(r)e^'(^(p,)) 
^ifM v'A2(r) - P(r) 

x^(r) =x^(0)-M^(j9,) r dnXin)+e^{u{ps)) f drAAn). (2.7) 

JO JO 

While the Dirac multipher A(r) is determined by the gauge fixing — r ^ 0, the 3 
Dirac's muhiphers A(r) describe the classical zitterbewegung of the origin x^(r) = 2;^(r, 0) 
of the coordinates on the Wigner hyperplane. Each gauge-fixing x{'^) ~ to the 3 first class 
constraints ^ (defining the internal rest-frame) gives a different determination of the 
multipliers A(r). Therefore it identifies a different worldline for the covariant non-canonical 
origin x^P'^{t) which carries with itself the definition dcom of the internal 3-center of mass, 
conjugate with /5+ 

Let us remark that the constant a;^(0) [and, therefore, also x^(0)] is arbitrary, refiecting 
the arbitrariness in the absolute location of the origin of the internal coordinates on each 
hyperplane in Minkowski spacetime. The origin x^{t) corresponds to the unique special 



Obviously each choice x{t) leads to a different set of conjugate canonical relative variables. 



i=l 



N 

Msys] - A(r)^Ki, 
1=1 



i'sir) = 
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relativistic center-of- mass-like worldline of Refs. [|T^] which unifies previous proposals of 
Synge, M0ller and Pryce quoted in that paper. 

The only remaining canonical variables describing the Wigner hyperplane in the final 
Dirac brackets are the non-covariant canonical coordinate a;^(r) |^ and p^. The point with 
coordinates x^(r) is the decoupled canonical external ^-center of mass of the isolated system, 
which can be interpreted as a decoupled observer with his parametrized clock {point particle 
clock). Its velocity a;^(r) is parallel to p^, so that it has no classical zitterbewegung. 

The relation between x^(r) and a;^(r) {a is its 3-location on the Wigner hyperplane) is 



esip^s + ^s) 



(2.8) 



After the separation of the relativistic canonical non-covariant external 4-center of mass 
a;^(r), on the Wigner hyperplane the N particles are described by the 6N Wigner spin-1 
3- vectors 57i(r), Ki(r) restricted by the rest-frame condition k+ = YJi=i ~ 0. 

Inside the Wigner hyperplane, three degrees of freedom of the isolated system Q become 
gauge variables. To eliminate the three first class constraints k+ ^ the natural gauge 
fixing is x(t) = (^com = g+ ~ such that g+ ^ implies Af.(r) = 0: in this way the internal 
3-center of mass is located at the origin a;^(r) = z^{t, a = 0) of the Wigner hyperplane. 

The various spin tensors and vectors are 0] 



V 

S 1 



= K(Ps)6'^(n(p.)) - u''{ps)e>^{u{p,))]Sr + e^{u{p,))e''{u{p,))Sr ^ 

N 



+ 



e^HPs)K{ps) - e^iu{p,))u^ips)\ E^rV^'c^ + + 

i=l 

N 

e'r{u{ps)K{u{Ps)) - e';{u{psM{u{ps))\T.^l<, 



i=l 



sr = ^MPs)K{u{Ps))sr 



16c 



See Refs. |2C,2T| for the definition of this concept in general relativity. By using the interpreta- 



tion of Ref . [ 20 1 , also the special relativistic limit of the general relativistic Dixon centroid of Ref . 



1 18 1 gives the centroid x'^{t): it coincides with the special relativistic Dixon centroid of Ref. |17] 



defined by using the conserved energy momentum tensor, as we shall see in Section IV. 

^"^It describes a point living on the Wigner hyperplanes and has the covariance of the little group 
0(3) of timelike Poincare orbits, like the Newton- Wigner position operator. 

^^As already said, they describe an internal center-of-mass 3-variable acom defined inside the 
Wigner hyperplane and conjugate to k+; when the (?com are canonical variables they are denoted 
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N 



N 



1 



i=l 



i=l 



QIJ ^1''^ S^^ S^^ 



PsM'p'^s - s^'p^) + ^JpKsrp': - stpt> 

Ps + \f^s' 



N 



s = s = J2vi ^ f^i 

i+1 



N N-1 
^ ^fji X Ki - f]+ X = "Y PaX TCa 



(2.9) 



i=l 



a=l 



Let us remark that while L'^'^ = x^p^ — x'^p'^ and S^'^ are not constants of the motion 
due to classical zitterbewung, both L'^'^ = x'^p^ — x^Ps S^'^ are conserved. 

The canonical variables x^, for the external 4-center of mass, can be replaced by the 
canonical pairs ||22[|p^ 



Ps ' '^S Ps ' 



-. / ~ Ps_ ~o\ 



Ps 



Ps 



(2.10) 



with the inverse transformation 



x: = ^Y7pats + ^ 



Xs ~l" (Tg + )^S5 



p: = es^Jl + kl 

Ps = ^sks- 



(2.11) 



This non-point canonical transformation can be summarized as [e^ — Mgys ^ 0, k,^ 



Ef=i f^i ^ 0] 





Vi 


p'i 





f^s 


Zs 


Vi 


T 

s 


— # 

kg 





(2.12) 



The invariant mass Mgys of the system, which is also the internal energy of the isolated 
system, replaces the non-relativistic Hamiltonian Hrei for the relative degrees of freedom. 



19 



It makes explicit the interpretation as a point particle clock. 
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after the addition of the gauge-fixing — r ^ 0: this reminds of the frozen Hamilton- 
Jacobi theory, in which the time evolution can be reintroduced by using the energy generator 
of the Poincare group as Hamiltonian 0. 

After the gauge fixings Tg — r ~ 0, the final Hamiltonian and the embedding of the 
Wigner hyperplane into Minkowski spacetime are 

Hd = Msys - A(r) ■ K+, 
z^{t, a) = <(r) + e^{u{ps))a' = <(0) + u^{p.)t + e^^(«(p.))a^ 



with 



i's{r) = + K(r), Ho] = u^iPs) + e^{u{ps))Xr{r), (2.13) 

where a;^(0) is an arbitrary point. 

The particles' worldlines in Minkowski spacetime and the associated momenta are 

x^{t) = z'^ir, v,{t)) = x^{t) + e^.{u{ps)H{T), 

pfir) = ^nq + R'i{T)u^{p,) + e':{u{p,))Ki,{T) epl = ml (2.14) 
The external rest-frame instant form realization of the Poincare generators P] with non- 
fixed invariants ep] = ~ M'^y^^ —ep^S^ ~ —eM^y^S , is obtained from Eq.(|2] 



s ~ -^sPs •^st's ^ '^s 1 



p: = ^el + pl = e,^l + kl ^ ^M^y^ + pl = M,,, Vl + ^1 

N 

1=1 

I 1 ^ 

^s + ^ei + pi i=i 

1 + k^^zl '-^= ~ 5X - xUm^^+jP, ^% ' (2.15) 

^rn ^ " M,ys+,/M^y,+pi 



^'^It implies A(r) = — 1 and identifies the time parameter r with the Lorentz scalar time of the 
center of mass in the rest frame, Tg = ps ■ Xg/Mgys', Mgys generates the evolution in this time. 

^"^See Refs. |23] for a different derivation of this result. 

^■^As in every instant form of dynamics, there are four independent Hamiltonians p° and J°*, 
functions of the invariant mass Mgyg] we give also the expression in the basis T^, e^, Zg-, kg. 
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On the other hand the internal reahzation of the Poincare algebra is built inside the 
Wigner hyperplane by using the expression of S^^ given by Eq. (|2.9| ) ^ 



N 



Msys = Hm = ^\lml + K. 

i=l 

N 

K+ = j2^i{^ 0), 



N 



J = ^ r/^ X r = 3'' = le™''^™ = SI 



i=l 

N 



^ = - E v^' + = -MsysR+, = r'' = SI 



i=l 

n = Mf^, -kI^ > 0, 

2 „ ^2 



= -e(M4, - Ki)S, ^ -eM^y.S,. (2.16) 



The constraints es — Mgys ~ 0, k+ ~ mean: i) the constraint — Mgys ~ is the 
bridge which connects the external and internal realizations 0; ii) the constraints k+ ~ 0, 
together with imply a unfaithful internal realization, in which the only non-zero 
generators are the conserved energy and spin of an isolated system. 

The determination of q+ for the N particle system was done with the group theoretical 
methods of Ref. in Section III of Ref. |1]]. Given a realization of the ten Poincare 
generators on the phase space, one can build three 3-position variables in terms of them 
only. For N free scalar relativistic particles on the Wigner hyperplane with psys = ~ 
and by using the internal realization (|2.16| ) they are: 

i) a canonical internal center of mass (or center of spin) g+; 

ii) a non-canonical internal M0ller center of energy 

iii) a non-canonical internal Fokker-Pryce center of inertia 

It can be shown [|l| that, due to k+ ~ 0, they coincide: ^ _R+ ^ 

Therefore the gauge fixings x{t) = gV ~ R+ ^ ^ imply A(r) ^0 and force the 
three internal collective variables to coincide with the origin of the coordinates, which now 
becomes 



X 



i^-+)^(T,)=x^(0)+M^(p,)r,. (2.17) 



As shown in Section IV, the addition of the gauge fixings x(r) = g+ ^ i?+ ^+ implies 



that the Dixon center of mass of an extended object |T8i and the Pirani [25| and Tulczyjew 



^^This internal Poincare algebra realization should not be confused with the previous external one 
based on S^'^; Ii and W'^ are the two non- fixed invariants of this realization. 

^^The external spin coincides with the internal angular momentum due to Eqs.( AlT| ) of Ref. [l|. 

^^As we shall see in the next Section, A ~ is implied by the natural gauge fixing g+ ~ 0. 
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26| , |21|p7| centroids all simultaneously coincide with the origin x^(r). 



The ea;terna/ realization ( p.l5|) allows to build the analogous externa/ 3- variables g^, Rg, 
Ys- It is then shown in Ref. ^ how to build the associated external 4- variables and their 
location on the Wigner hyperplane 

Ps 

= (x- ^ + (T, + ^)ks) =x^ + e(l(«(p.))a", 



a 



Y 



= x^,+eMPs))a^, 
R^ = (x:; Rs) = 

' ''"e,M°(p,)[l+M°(p,)] 
= ■ Xs = u{ps) ■ Xs = u{ps) ■ Ys = u{ps) ■ Rs, 

erMPs))[MPs)S7 + sr] _ 



[1 + U-{ps)] 



orr I Ps ^ r>r , (Ps) 

-^s H rr- — — — = ^3^4- + 



es[l + u"iPs)] l + u-ips) 

r S:-'u%Ps) Sl'u^Ps) 

'^^+ + i + uo{p,) ^ i + uo{p,y 



ts[l + U"{ps)] 



1 + U°{ps) 



^^See Ref. for the application of these methods to find the center of mass of a configuration 
of the Klein-Gordon field after the preliminary work of Ref. |29| on the center of phase for a real 
Klein-Gordon field. 
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esu''{ps)[l+u''{ps)\ 

^~r_ S^U^p,) ^^^r ^ [l-U%Ps)]SYu^{p,) ^ 

U°{j),){l+U°{j),)] " + U°{ps)[l + U''(j)s)\ 

_ [l-U°{ps)]S:'u'{Ps) 
U''{ps)[l+U^{Ps)] ' 

^ ^ y^/.^ (2.18) 

namely the external Fokker-Pryce non-canonical center of inertia coincides with the origin 
xf+^^^ir) carrying the internal center of mass. 
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III. THE ENERGY-MOMENTUM TENSOR ON THE WIGNER HYPERPLANE 
AND DIXON'S RELATIVISTIC MULTIPOLES. 



A. The Euler-Lagrange Equations and the Energy-Momentum Tensor of 
Parametrized Minkowski Theories. 

The Euler-Lagrange equations associated with the Lagrangian are (the symbol '=' 
means evaluated on the solutions of the equations of motion) 

dL _^dL _ 

_rl^ll ^^^^ f) 9rr + grsVi , ^ x 

where we have introduced the energy-momentum tensor [here rjfir) = (l;r/j(r))] 



Due to the delta functions the Euler-Lagrange equations for the fields 2;'^(r, a) are trivial 
(0 = 0) everywhere except at the positions of the particles. They may be rewritten in a form 
valid for every isolated system 

c^^T^^^ ^ - ^9^[^4]r^^. (3.3) 
V 9 

When dAly/gzj^] = as it happens on the Wigner hyperplanes in the gauge g+ ~ and 
T, — r ^ 0, we get the conservation of the energy-momentum tensor T^^, i.e. OaT^^ = 0. 
Otherwise there is a compensation coming from the dynamics of the surface. 
On the Wigner hyperplane, where we have 

xr(r) = z>^{T,Ur))+z^{r,m)vl{r) = i;^(r) + (n(p,))r)r(r), 

= 9rT{T, Vi{r)) + 2grr{T, Vt{T))fi{T) + Qrsir, Vi{^))ViX^)vtiT) = 



= y^^JT^K^M'^(p,) + e^(M(p,))<(r) ^ p2 _ 

N 

p^= d'ap,{T,a)^Y.P'i^)^ (3-4) 

1=1 

the energy-momentum tensor T^^{t, a) has the form 
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-1 \lil{r) + 2xs,{T)^r{u{ps)) - fkir) 

^ mi?)[(r) 



i=i ^Jxl{T) + 2i,^(r)6(^(w(p,)) - r^(T) 

T^^lr, a) = - E - m) I ^-^'^(")^-^(^) . (3.5) 

-1 V^2(^) + 2i;,^(r)e^f(«(p,)) - 7^(r) 



B. The Energy-Momentum Tensor of the Standard Lorentz-Covariant Theory. 

The same form is obtained from the restriction to positive energies "^"^ of the energy 
momentum tensor of the standard manifestly Lorentz covariant theory with Lagrangian 
Ss = I drLsir) = — J2i^=i I dT^Jxlir). On the Wigner hyperplanes with — r ~ we 
will get 



i:r(ri)i;r(ri) 



= E^^ / rfr,^iA^f^5^(x,(rO - z{r,a)) = 

i=l \Jxi[Ti) 

^e^Au{ps))eUu{PsW''{T,a). (3.6) 
On the other hand, from the restriction of the standard theory we get 

1) On arbitrary spacelike hypersurfaces 



1=1 



N 



i=l 

[z>^{tu ffiin)) + <(ri, ffiiriMin)] 

[<(ri, rfi(ri)) + <(ti, f]i{ri))r]l (ti)] = 

N r dr 

/ r^7^ ^X^(^i^ Vi{n)) - z{t, a 

i=l -I Jxi[Ti) 



Pi = mi 



[z>:{nMri)) z^inMri)) + 

+ (^,^(n,7f,(ri))<(Ti,rf,(Ti)) + <(Ti,r/,(Ti))<(Ti,ry;(Ti)))77[(Ti) + 
- > 0, e(EiIiPn' = <T.lim^f > 0. 
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+ < (n , ffi (ri ))z''^ (n , ff, (n) )?)[ {n)fit (n) 

N 



+ <(r, r/,(r));.,^(r, Ur)M{r) + (r, r/,(r))<(r, r7,(r))7)[(r)7),^(r) 



[4(r, r7,(ri)) z,^(r, r/,(r)) + (^(r, r/,(r))<(r, r/,(r)) + 
+ <(r, r/,(r))<(r, r7,(r)))7)[(r) + <(r, r/,(r))<(r, r7,(r))7)[(r)7)Kr) 



since det \z'X\ = ^/g = Ji{9tt - Y^arrgrs), 7 = \detgr 



2) On arbitrary spacelike hyperplanes, where it holds 
^^(r, a) = x^{t) + b^^{T)a\ xf (r) = x^{t) + b^^rK^r), 



z^{t, a) = Kir), z>^{t, a) = x^(r) + 6^:(r)a" = F/^ - g^^z^^, 



9tt 
9rs 

9 = 9tt + J29 



7 = 1, 



2 

rr' 



mm2 



g^^ = + b^^anV, 9^' = 9''9rr = br,{x^, + &!la")/[/^(a:^ + 



we get 



«m2 



1^, 



=1 V^(^, ff)y9rTir, ff) + 2g^rir, a)r]i{T) - t]^ (r) 
(x^(r) + KirMirmx:iT) + b:{rK{r)) + 

+ KirK{r)rj:iT)fjt{T) . 



3) On Wigner's hyperplanes, where it holds 
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9rr = is^Cr 9 



g = [Xs{t) ■ U{ps)]'^, Qrr = X^, 

g^^ = l/[xsXiPs)]\ 9^' = Xs,e^.Hps))/[xsXiPs)]^ 
grs = _srs ^ x,^e^{u{ps))x,,e''^{u{p,))/[xs^,u''{ps)]\ 



rs ^rs} 



ds^ = xl{T)dT'^ + 2is(T) ■ er{u{ps))drda'^ — da'^ 



(3.10) 



we get 

Tn^'s{r) + e'rHPs)y] 



N 

E 



rrii 



6'{a-v.{T))[x^iT)x:iT) + 

+ {x^{T)e-:.{u{Ps)) + i:{r)e^HPs)))vKr) + 
+ <(«(P.))<(«(P.))^[(r)^-(r) 
1 



N 

E 



TTli 



isir) ■ u{ps) v'i;2(^) + 2xsp{T)e?{u{p.))rjl{T) - f^(r) 
x^(r)xr(r)53(a-r/,(r)) + 

+ (i^M<(«(Ps)) +xr(r)6tf(«(p.)))53(a-7^.(r))^[W + 

+ e^{u{p,)YXu{Psy)^\a - r/,(r))7)[(r)7),f (r)] . (3.11) 

Since the volume element on the Wigner hyperplane is u'^{ps)d^a, we obtain the following 
total 4-momentum and total mass of the N free particle system (Eqs. ( |2.7] ) are used) 

P^= f d'aT^''[x^,{r) + e^MPs))cr'']uM = 
rrii 



N 

Sv^)-K(r) + A(r)P 
[-X{rmPs) + mr) + y{r)]e'^HPs)) 

N AT 

Ur)=-i E ymh' + nKT)u^{Ps) + <(r)elf(n(p,))l = (^) = P^, 



i=l 



i=l 



N 



M,y, = P^u^{p,) = -Xir)Y^ 



rrii 



.=1 VA2(r)-[r/,(r)+A(r) 



N 



lA(r)=-i Y.y^f + ^ii'r 



(3.12) 



i=l 



which turn out to be in the correct form only if A(r) = —1. This shows that the agreement 
with parametrized Minkowski theories on arbitrary spacelike hypersurfaces is obtained only 
on Wigner hyperplanes in the gauge Ts — t ^ 0, which indeed implies A(r) = —1. 



21 



C. The Phase-Space Version of the Standard Energy-Momentum Tensor. 



The same result may be obtained by first reformulating the standard energy-momentum 
in phase space and then by imposing the restriction (r) = 

^mf - 7^«(r, r/i(r))Ki^(r)Ki,(r)P(r, r/i(r)) - /«i^(r)7''"(r, r/i(r))z^(r, r/i(r)) 



mf + /??(r)M^(ps) + i^i^ri'^iPs)) ^^[the last equality refers to the Wigner hyperplane, see 
Eq.(^]: 



i=i mi^g{T, a) 
J| on Wigner' s hyperplanes 



N 



T^'^[^fW + 6f(«(p.)K] = E 



i=l 
N 



x l + 2xsi3efu{u{ps))r]i - rj 
mi^Xs ■ u{j)s) 



i=l 



nii^Xs ■ u{ps) 
[{m^ + k^{t))u^{Ps)u''{Ps) + 
+ kl{r)^m^ + Kj{r){u^{ps)e^MPs)) + u\ps)e^Au{p,))) + 
+ <{r)Kl{T)e>:{u{ps))e:{u{ps))\ = 



N 



E'^'('?-^iW) 



\\T)-%{r)+\{r)f 



i=l 



nii 



mf + Kf{T)u^\ps)u''{ps) + 



+ kl{T)[u^{ps)e^MPs))+u''{ps)e^{u{ps))} + 
<(r)<(r) 



+ 



mf + Kf(T 



-.e'f{u{p,))e:{u{p,)) 



(3.13) 



Now the total 4-momentum and total mass are 



N ™2 



E 

i=l 



xl + 2xsf,e^MPs))v^-%, Jmf + Kf{T) 



Xs ■ U{ps) 



mi 



^^In this way we are sure to have imposed the restriction to positive energy particles and to have 
excluded the other 2^ — 1 branches of the total mass spectrum. 
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[sjml + ^,{t)u^{Ps) + n\{T)e^Au{Ps))\ 



TV 



M,,, = P^uM = E ' ^ ' j V ^f + (t). (3.14) 

These equations show that the total 4-momentum evaluated from the energy-momentum 
tensor of the standard theory restricted to positive energy particles is consistent with the 
description on the Wigner hyperplanes with its gauge freedom A(r), A(r), only by working 
with the Dirac brackets of the gauge Tg = t, where one has A(t) = —1 and 

x^g{T)=x>:{0)+TU^{Ps) + e'^{u{ps)) rdnXrin), (3.15) 

Jo 



because rrii/yl — [f/j(r) + A(r)]2 = ym-c^ + Kf(r). 
Therefore, for every A(r), we get 



T^''[x^g{Tg) + 6^(«(p.))a1 = e^^{u{pg))eUu{Ps))T^''{Ts, a) = 

= E ^'i^ - ViiTs)) ym^ + K^iTg)u^^iPsKips) + 



i=l 



+ A;[(T,)(n^(p.)<(«(P.)) +«''(P.)e(^(«(p.))) + 



N 



T^^iTs: ^) = E ^'(^ - UTs))^Jml + (T,)), 

i=l 

AT 

r-(T„a) = E^'(^-^^(7^3))<(r.), 

= = M«^(p,) + e':{u{ps))K\ ^ Mu^'ip,), 

N 



i=l 



T>^Vs{Ts) + e^MPs)K] MPs) = e'X{u{p,))T^^T„ a) = 



29Namely we get = pi^ and Msys = Eili \frr^^r^Jj)- 
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N 



[^mjc^ + ^,{T,)u^{ps) + ^\{Ts)e>^Au{Ps))l 



N 



1=1 



(3.16) 
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IV. DIXON'S MULTIPOLES ON THE WIGNER HYPERPLANE. 

In this Section we shall define the special relativistic Dixon multipoles on the Wigner 
hyperplane for the N-body problem^°. We list the non-relativistic multipoles for N free 
particles in Appendix A by comparison. 

Let us now consider an arbitrary timelike worldline w^{t) = z'^{t, 
e^(M(Ps))?7''(r) and let us evaluate the Dixon multipoles Jl^ on the Wigner hyperplanes in 
the natural gauge with respect to the given worldline. A generic point will be parametrized 
as 

(4.1) 

so that SZfj_{T, (y)u'^{ps) = 0. 

For 57(r) = we will get the multipoles with respect to the origin x^^+^^(r). 

A. Dixon's Multipoles. 

Lorentz covariant Dixon's multipoles and their Wigner covariant counterparts on the 
Wigner hyperplanes are defined as 

= e^;(n(p.))...e^f:(«(p.))e!;(n(p.))e^(«(p.))g^.--^^(T„rf) = 

= e'X{uips))eUuiPs)) I d^aSz'^^{Ts,a)....6z^'^{T,,a)T^''{Ts,a) = 
= e'^l{uips))...e^.:{u{p,)) 

r ^ I 

u'^iPs)u''iPs) EivViTs) - v''{Ts)]...[vl"iTs) - r/-(T,)]v/m| + K^iT,) + 
+ e^{u{ps))e^MPs)) 

i=i ^m^ + i^XTs) 
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Previous studies of multipole theories of point particles can be found in Refs. 
^,E7l^j33,|3l|,l33|-|35i. 



^^See this paper for the previous definitions given by Bielecki, Mathisson and Weyssenhof 
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N 



1=1 



N 



N 



N 



= J d'a5z'''{T,a)...Sz''-{T,a)T%[x'f+^''{T,) + e'^{u{p,))a^] 
^e':i{u{p,))...eMu{p,)) 



N 



i=i sjmi + Ki{Ts) 

= 6(f;(u(p,))...6$^;(u(p,)) 

i=l 

Related multipoles are 

= e^;(«(p,))....6^;(«(p,)) 



i=l 

+ KKr)u^iPs) + kliT,)e>^iuip,))] , 

u,Mpi}.'-^-^iT„ff)=0, 

The inverse formulas, giving the multipolar expansion, are 



26 



oo ri...r„AB/^ -.x 



n=0 ^• 



n=0 



-6\a-rfiT,)). (4.4) 



The quantities g^,i-""""(T„ r/), g;^-""'^"(T„ rf) = g5,i-''""'^(T„ r/), g;^-""™(T„ rf) are the 
mass density, stress tensor and momentum density multipoles with respect to the worldhne 
w'^lTg) (barycentric for fj = 0) respectively. 



B. Monopoles. 

The monopoles correspond to = '^^ and have the following expression 

i=i Jmf + Rf 



TV 



1=1 V de 



^^They are ff independent; see Appendix C of Ref. and Appendix A for the non-relativistic 
Hmit. 



^■^In this Section we use results from Section V of Ref. [0, where the rest-frame (k+ = = 0) 
canonical relative variables with respect to the internal 3-center of mass were found by means 
of a classical Gartenhaus-Schwartz transformation [14|. This transformation is a sequence of 
canonical transformations depending on a parameter a. The final rest-frame relative variables 
are obtained in the limit a ^ oo. In Eqs.(|45|), (|]9|), ( |4.11| ), ( |4.13| ), ( |4.14| ) we use the fol- 



lowing symbols from Ref. |l|: Ki(oo) 



NEa=lla^^c,a, H 



M2 



Hi 



mf + Kf 



lailbiT^qa ' T^qb- Also the uou-relativistic limits, c — > oo, are shown. 
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N v^l..7V-l =! =! 

-\- 0^0^ i\ . — 

i=l J mi + N Y.de IdileiT^qd ' T^qe 



N 1..7V-1 AT 

QT^'iTs, ff) ^c^oo ^i^^ + 9 H H ^^^ga " T^qb + 0{l/c) = 

i=l ^ ah 1=1 

N 

= "^iC^ + Hrel,nr + 0(l/c), 

j=l 

q^{Ts,f]) = k'^ ^ 0, rest — frame condition {also at the non — relativistic level), 

l..N~l N ^ 

g^nT.,rf)^™ Y E^^^« + 0(l/c) = 

ab i=l * 
I..N-1 1..N-1 

= E ^aY« + o(i/c)= E w>;: + o(i/c), 



N 2 



N 2 N 1 

mi -A mi 



Q— >00 

1 = 1 

N 1 l..N~l N 

"^iC ~ 9 TTga ■ TTqb + 0(1/C) = 

i=l ^ ab j=l 

TV 

= Y^^C^ - Hrel,nr + Oil/c). (4.5) 

j=l 

where we used Eqs. (5.10), (5.11) of Ref. |l|] to obtain their expression in terms of the 
internal relative variables. 

Therefore, in the rest-frame instant form, the mass monopole is the invariant mass M = 



J2iLi Y + while the momentum monopole vanishes. 

C. Dipoles. 



The dipoles correspond to n = 1 



35. 



qlr^^iT,, fj) = 5^5?M[R\{T,) - r^^(T.)] + 5„X[E ^==(^^) " ^^m)?T ^1 

i=i ^mf + Ki 

N 



+ 



35t^i^^(T„,?) = 6Af^^(^x(p,))e:^(^z(p,))e^(n(p,))4i^^(T„r?), if^\T,,fj) = e^/ (n(p,))4^^^(r„ r/), 
ui/ifrp -f. _ ^ui/„,/^ w^iAt/ 
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i=i ^mf + Kf 

The vanishing of the mass dipole identifies the internal M0ller center of energy i?+ ^ 
g+ ~ y+ and therefore the rest-frame internal center of mass (f+: 

q^T^^iT^ff) = e;\{u{ps)W{Ts,ff) = 0, ^ r/(T,) = for a;^ = x^)^]. (4.7) 

The time derivative of the mass dipole identifies the center- of-mass momentum-velocity 
relation for the system when ff = 

^^^^^ = Kl- Mf{T,) -.,^^0 0. (4.8) 

The expression of the dipoles in terms of the internal relative variables when fj = = 
g+ = is obtained by using Eqs. (5.10), (5.11), (5.22), (5.25) of Ref. [H 

qTiTs^R+) = 0, 

TV N-1 

i=l a=l 
N-l 



E r u 
rqa qa 

a=l 
N-l 



a^oo / , rqa"qa 
a=l 
N-l 

'C— +00 

a=l 



a— >oo 



i=i Hi Hi 

^EE7.p:;^ + W-i?;)E^ 

i=l a=l i=l 
2^ 2^ [lai-laj)pq 



E (cV^E(7a.-7a,) 



- >Va] + N Zde Idjlej^qd " ^qe 



a=l ij J mf + N J2de IdileiT^qd ' T^qe 



9C \ r 



Efc=l 'Jmi + NY.de IdklekT^qd ' T^qe 



1..N N-l _ ]\T 1..N-1 

c^oo 



1 l..A^-l ^, ^, ^, Y^A^ 



abc i=l '''' 



N 2 1 l..AfAf-l o- 

q^/A{Ts,R^) = EK-^;)^ = 47EE(7a.-7a,)p;f^ 
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N-l 1..N 



a=l 



m] + N J2de IdjlejT^qd ■ 7f, 



qe 



mf + N Y.de IdileiT^qd ' 



X 



qe 



rrii 



Pa 



Ek=l yml + N J2de IdklekTTqd " T^qe ' 
\..NN-\ ^.jyL.Ar-l 
^c^oo XI H ^{lai - laj)p'a— IhilciT^qb " T^qc + 0(l/c) 



ij a=l lio^iiu 

1 N—l N v^Af 



abc i=l 



m 



Pln^b-^qc + 0{l/c). (4.9) 



The antisymmetric part of the related dipole pj.^^(Ts, 17) identifies the spin tensor. Indeed, 
the spin dipole Q is 



+ 



N 



+ YivKTs) - r/^(T,)]<(T,) e(f(«(p,))e:(«(p,)) - e^AuiPs))e^sHPs)) 



i=l 



5r(T,)[r/=0] = 5r = 



-1 Vl-r/,(T,; 



e^?(«(p.))er(«(p.))-e,^(«(p.))e^(«(p.)) 



1 

= -e';{u{ps))M[Rl{T,)-r^^^{Ts)] = -eUu{psWT^{Ts,v), 



u^{ps)Si^''{T, 



r/= 0. 



(4.10) 



This explains why ^(^(p ){Ts,f]) is also called the mass dipole moment. 
Therefore, x[^+'^^{Ts) is also simultaneously the Tulczyjew centroid p6| , pT| , p7| Qand, due 
to x^f+^'^{Ts) = u^{Ps)-i also the Pirani centroid |^ Q. Usually, in absence of a relation 
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See Ref. |17l for the ambigous definitions of p^,, Sjf given by Papapetrou, Urich and Papapetrou, 



and by B.Tulczyjev and W.Tulczyjev [38|-40|. 



^^Defined by Sjf {Ts)[Tf\Ui,{ps) = 0, namely with 5'°'' = in the momentum rest frame. 
^^Defined by Sl^ {Ts){f]\xsu = 0, namely with S""^ = in the instantaneous velocity rest frame. 
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between 4-momentum and 4-velocity they are different centroids 0. 

Let us remark that non-covariant centroids could also be connected with the non- 
covariant external center of mass and the non-covariant external M0ller center of energy. 



D. Quadrupoles. 

The quadrupoles correspond to n = 2 PI: 



N 

4 = 1 

i=i ^mf + Rj 

N 



i=l 

N 

ql?^-^-{Ts, R+) = J2iv? - R7)iv? - R7)H^ = 

i=l 

^ 1..N1..N-1 H H H 

ijk ab ^ 

N 

1..N1..N-1 , K^HjH , 
ijk ab M 



i=l 



J^Y. hai-laj){lbi-lbk)PaP?^J^-^, 

i=i Hi 



1 I..NI..N-1 K^K^HH 

= nYJ: (7a. - 7a,)(7M - lbk)p7p? '^^r . (4.11) 

ijk ab * M 

Dixon's definition of barycentric tensor of inertia follows the non-relativistic pattern 
starting from the mass quadrupole 

N 



g^,^^^--(T„r/) = YM'V^y^ + ^']iTs) + M[^^W'' - V^'R^ ' V^'^R^W, (4.12) 



1=1 



^^For instance, the so-called background Corinaldesi-Papapetrou centroid |41] is defined by the 
condition Sjf{Ts)[ff\Vi, = 0, where is a given fixed unit 4-vector. 

X'^"^' = e^f^i(n(p,))e^;2(^.(p,))e^(n(p,))e^(n(p,))g;i^'^^^, 

= e^^l{u{p,))e^^l{u{ps))e'X{u{ps))q^'''^^A, t^^^^ = (^z(p,))e^f^2 (^^(p,))^^^^^-, p>^'^^^ = 
e'^^{u{ps))ef^'^{u{p,M{u{ps))qJr''''''^. 
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u 

N 



=] 
1 



1..N 1..N-1 



ijk ah 

l..Af-l 1..N 

J2 {j;^Y. (7ai - 7ai) (7M - Ibk) 
ab ijk 



/m^ + ^ 'JdjlejT^gd ■ T^qe ml + ^ H IdklekT^qd " TTge 



de V de 

ri 



afe ijk 

^l..A^-l , -2 -2 Arv^l..Ar-l 



IcildiT^qc- T^qd ^ N Ecd 7ej7d?-7rgc " T^gd _^ 

c ^ 2m? 2m? 

, ^Ecd^~Scfc7dfc7rgc • TTgrf 1 ^ A^Eed^~Scfe7dfo7rgc • TTgrf X p/-| /2V 



l..Af-l 



r''^[qnr]+0{l/c). (4.13) 



In the non-relativistic limit we recover the tensor of inertia of Eqs. (|All|) . 
On the other hand, Thome's definition of barycentric tensor of inertia |H2 



IS 



i=i \ mf + K.f 



^ 1..N1..N-1 

]^ E E (7a* - laj){lbi - Ibk) 
ijk ab 

[p.-pV^^^-pM-"^^''^-^^^''^^^^ 



Hi{oc)U 



1..N-1 ^ 1..N 



E (]y E (7ai - laj){lbi - Ibk) 
ab ijk 
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mfJmj + N Y.de IdjlejT^qd ' T^qeJ mi + N Y.de IdklekT^qd " 



mf + N J2de IdileiT^qd " T^qe{Eh=l \j + ^ Ede IdhlehT^qd " T^qeY 

ri r2-| 
qarqbi 

{lai - laj){lbi - IbkjPqa 



[Pqa ■ PqbS'^^' ' p'^J.' 



\..N-\ \..N 

22 E 



+ 



ab ijk 

i + i( 



Nm? 



PqaPqbl ^ 



1..N-1 
cd 



IcildiT^qc ■ T^qd . l^cd IcjldjT^qc ' T^qd 



2mf 

IckldkT^qc ■ T^qd 



+ 



1 AiVE^/-^ 



2ml 



m 



E 

h=l 



2m'j 

IchldhT^qc • T^qd 



+ 



rrih 



+ 0(l/c2 



l..Af-l 



'[g;,] + 0(l/c). 



(4.14) 



In this case too we recover the tensor of inertia of Eg. ( |A11|) . The Dixon and Thorne 
bary centric tensors of inertia differ at the post-Newtonian leveh 



Trir2 
dixon 



jrir2 (rn \ 
^thorneK-'' s) 



1 ^-^-^^-^ m-m^ 

- E 2^J^ila^-laJ){r^bi-lbk) 
ab ijk 



Pqa ■ PqbS'""'- 



PqaPqb 



l^cd IcildjT^qc ■ T^qd 

rrii 



0(l/c^ 



E. The Multipolar Expansion. 

In the next Section further types of Dixon's multipoles are analyzed. It turns out that 
the multipolar expansion (|4.4|) can be rearranged with the help of the fact that Hamilton 



equations imply d^T^^ = so that from Eqs.( ^.8|) we get the multipolar expansion 
T'^'^lxf+^^iTs) + e^.{u{ps)y] = T^-'lwfiTs) + e^(n(p,))K " V^iT^))] = 

= u^'^iPsyJiuipMS^M + 6^Kl]5%a - r/(T,)) + 

oo (T\n an 

+ E ^/^-^"^^(T., r/)6- - m)), (4.15) 

where for n > 2 and r/ = /^^-^"^"(T,) = 4a*i-/^"-iIa'I/^")-(t,), with J^'-^-^''P''{Ts) 

being the Dixon 2'^^"' -pole inertial moment tensors given in Eqs.( |5.10D . 

The equations d^T^'^ = imply the Papapetrou-Dixon-Souriau equations of motion 
T8|^ , |4^ for the total momentum Pj^iTg) = e'X{u{ps))qT^ {Tg) ^ and the spin ten- 



sor Sj^{Ts)[fj = 0] restricted to positive energy particles [see Eqs. (|5.4| ) and (|5l7| )] 
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^ 2P^^(T,)u%,) = 2Klel^{u{p,))u%,) ^ 0. (4.16) 
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V. MORE ON DIXON'S MULTIPOLES. 



In this Section we shall consider multipoles with respect to the origin, i.e. with rf = 
[we use the notation t^^-^"'^"(T„ 0) = t^^-'^"'^"(T,)] 

As shown in Ref. Il^, if a field has a compact support W on the Wigner hyperplanes T,wt 
and if f{x) is a C°° complex- valued scalar function on Minkowski spacetime with compact 
support 0, we have 

<T^"',f> = I <fxT^"'{x)f{x) = 



oo (_j)n 
n=0 ^• 

r r d'^k °° (—iV^ 



(2^) 

and, but only for f{x) analytic in W Ijl^ Q, we get 

r °° 1 (9"ffr 



\x=XsiTs), 



T-(x) = E ^ ^/ / dTs5\x - x.(T.))t--— (T.). (5.2) 



For a N particle system this equation may be rewritten as Eq. 
For non-analytic functions f{x) we have 



/• r d^h ~ °° ( — i']"- 



n=N+l 



and, as shown in Ref. |17], from the knowledge of the moments tj^"''^"'^[Ts) for all n > we 
can get T^^{x) and, therefore, all the moments with n < N . 
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So that its Fourier transform f{k) = J d'^xf{x)e'^^'^ is a slowly increasing entire analytic function 



on Minkowski spacetime (|(x° + iy' 



o\qo 



Ax 



< C pa-o\y°\- 



■+a.3\y^\ 



positive integers for every /i and Cq^,„q^ > 0), whose inverse is f{x) = J ^f^/(^)e 



an > 0, 
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See this paper for related results of Mathisson and Tulczyjev |36,30,^,2f 
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The Hamilton equations imply imply (we omit the dependence on 57 = of the multi- 
poles) 



If we define for n > 1 [Si^ = 2pJf = 2cj^^ 



0, forn = 0, 

-nu^^'{ps)pii?-^"^^{Ts) + nti[^'-^"^^{Ts), n>l. 



(5.4) 



= e(f(«(p.))....e:f:(^(p.))e^(«(p.))g---^^(T.), 



1 



{u{p,))....e^;MPsm'-'"'iTs) = ^«^(Ps)??-'^"^^(T.) + e^f(«(p.))gf-^"'-)^(T,) 

n + 1 

= [e$f^Htz(p,))..<"e:;(«(p,)) - 
-6(f(n(p.))...6^:(«(p.))e^(n(p.))]g5.--^^(T,), 

J^^l-l...^ln^l) , 



'Ml 



T) 

+ 1 

+ 6^f(n(p,))[g-----(T,) - gj:— )-(T,)], 



(5.5) 



and then for n> 2 



{^ll...^^n)(^lu) 

T 



(T,) = t^^-'^"^"(T,) - 



^ T ST S 



+ 



+ 



(T.) = 

n + 1 



n 



{lJ.l...flrifJ.)u 



(Ts) = 0, 



^^In Ref. [0 this is a consequence of d^T^^ = 0. 
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T) — 1 

n 

[(n - l)g?,-^"^^(T,) + gf-^'"^)^(T,)] + 



+ e^Au{ps))e:MPs)WT (Ts)- 



n 

Eqs.(O) may be rewritten in the form 
1) n = 1 

2 dig 

dTg 



=1 \/mf + /tf 



5^^(T,) = 2p^^{Ts)u^\ps) = 2k14^{u{ps))u''\ps) ^ 0, 



dT. 



2) n = 2 [^c/ent^ty = 4'''')'^ + 4""^^ + 4'^'^)''] 



3) n > 3 



This allows to rewrite < T^'^, / > in the following form |T7| 



oo ^ _■ 



n=2 ^ 



37 



with 



TfJ.l...fJ,nfJ.l' 



2 



n — \ 



n 



n + 1 



_|_ pP-nf^) ^i^ I _|_ 



n 



ri..r„AB 
HT 



(Ts) 



ri •••"=r„_i"=r„'^A '^ri •••'^r„-i^rn'^A 



te: 

n 



(f;..<"6!f-e(^^..e^:e!r)n'^)(p,)l(n(p,)) 



r(/xi.../i„/i)i/ 



(T.) = 0, 



Ml 



n + 
n 



(5.9) 



For a N particle system Eq.(5^) implies Eq. (f4.15 ). 
Finally, a set of multipoles equivalent to the j^i--^"^'^ jg 44 



for n > 



^^ll...^ln^lupcr 



1 



r/il...At„[At[pH'^] 



^|J^...^lr^[^l[pv]a] ^JT^ _ 



n + 1 



u 



T 



44^[/x[p!^]o-] "^j/ l(^j^pi^p<J _ j^uppcF _ j^pcrup _j_ j^uupp-j 
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n 



+ u^%Ps)e:i{u{pM'HPs))e'i{u{ps)) 



(n + 4) {3n + 5) linearly independent components, 



u^,{Ps)J^'-''"''""'{Ts) = jA^i-Mn-i{/.nM-)p<x^y^^ _ forn> 1, 



A(n - 1] 



n + 1 



-J. 



{fj,l...fl„-l\fl\fln)u 



T 



for n > 2, 



e-(«(p,))....6-(n(p,))J^^ 



...^in^ivpa 



1 



n + l 

+ u^'{Ps)e:\u{pM'HPs)VMPs)) 



(5.10) 



The j^i--/^"/^'^''"' are the Dixon 2^^'^-pole inertial moment tensors of the extended system: 
they (or equivalently the jw - Mn/i^^'g^ determine its energy-momentum tensor together with 
the monopole Pj< and the spin dipole Sq^ . The equations dfj_T^^'^ = are satisfied due to the 
equations of motion ( |5.7|) for and 5*^*^ without any need of the equations of motion for 



the J^^^■■^^"^^^P^ _ When all the multipoles ^re zero (or negligible) one speaks of a 

pole- dipole system. 

On the Wigner hyperplane the content of these 2'^^'^ -pole inertial moment tensors is re- 
placed by the Euclidean Cartesian tensors q^'"'^"'^'^ , q^'"^"^^, q^'"^"^"^. As shown in Appendix 
B we can decompose these Cartesian tensors in their irreducible STF (symmetric trace-free) 
parts (the STF tensors). 

The multipolar expansion (|4.15| ) may be rewritten as 



T'^'^ixf )'3(r,) + (w(p,))aT = T^-'lw^iT,) + («(p,))(a'- - v''{T,))] 



= uHPsylHp^S^M + 5^Kl]5%a - rf(T,)) + 
+ ls'T^'{TsMu^\ps)e;{u{ps))^6'{a - r/(T.)) + 



l)"rn + 3 



n\ 



n + l 



u^^{Ps)u^{Ps)qT'-'"^^{T,,fj) + 



+ E- 

ra=2 

+ -K(p.)e,^(n(p.)) + n-(p.)6(f(n(p.))]g^,-^"'^-(T„r/) + 
n 
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The so called Papapetrou- Dixon- Souriau equations given in 
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e'sMPs)KMPs))W---''-''''{Ts,v^ 



n 



(5.11) 
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VI. CONCLUSIONS. 



In this paper we have completed our study of the relativistic kinematics of the system 
of N scalar positive-energy particles in the rest-frame instant form of dynamics on Wigner 
hyperplanes orthogonal to the system total 4-momentum initiated in Ref. IQ. 

We have evaluated the energy momentum tensor of the system on the Wigner hyperplane 
and then determined Dixon's multipoles for the N-body problem with respect to the internal 
3-center of mass located at the origin of the Wigner hyperplane 0. In the rest-frame in- 
stant form of dynamics these multipoles are Cartesian (Wigner- covariant) Euclidean tensors. 
While the study of the monopole and dipole moments in the rest frame gives informations 
on the mass, the spin and the internal 3-center of mass, the quadrupole moment give the 
only (even if not unique) way to introduce the concept of barycentric tensor of inertia for 
extended systems in special relativity. 

Finally, let us observe that, by exploiting the canonical spin bases of Refs. []TB|,P^, after 
the elimination of the internal 3-center of mass {q^ = = 0) the Cartesian multipoles 
g^^-''"^^ can be expressed in terms of 6 orientational variables (the spin vector and the 
three Euler angles identifying the dynamical body frame) and of 6A^ — 6 {rotational scalar) 
shape variables. 



So that the origin is also the Fokker-Pryce center of inertia and Pirani and Tulczyjew centroids 
of the system. 
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APPENDIX A: NON-RELATIVISTIC MULTIPOLAR EXPANSIONS FOR N 

FREE PARTICLES. 



In the review paper in Ref. [0 there is a study of the Newtonian multipolar expansions 
for a continuum isentropic distribution of matter characterized by a mass density p{t,a), 
a velocity field U^(t,a) Q and a stress tensor a^'^(t,(T). If the system is isolated, the only 
dynanical equations are the mass conservation and the continuum equations of motion re- 
spectively 



dp{t,a) dp{t,a)U^{t,a) _ 

dt da^ 
dp{t,a)U^{t,a) djpU^U' -a'']{t,a) ^ 

di da' ~ 



0. (Al) 



We can adapt this description to an isolated system of N particles in the following way. 
The mass density 

N 

p{t,a)=Y.m,5\a-f],{t)), (A2) 

i=l 



satisfies 



dt ^ ^ " da'' 



while the momentum density 



N 



p{t,a)U'{t,a) = J2mifJ,{t)6%a -ffiit)), (A4) 



i=l 



The associated constant of motion is the total mass m = J^iLi- 

If we define a function ({a,f]i) concentrated in the N points ffi, i=l,..,N, such that 
C{a, ffi) = for (J 7^ rfi and ({ffi, fjj) = 6ij then the velocity field associated to N particles 
becomes 



AT 



U{t,a)=Y.mC{ff.m)- (A5) 
1=1 

The continuum equations of motion are replaced by 



'^'^p{t,(j)U{t,a) is the momentum density. 

^^This can be taken as the definitory equation for the velocity field, even if strictly speaking we 
do not need it in what follows. 

^^It is a limiting concept deriving from the characteristic function of a manifold. 
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Pj ft ^ ^ 

def d\pU^U^-a^^]{t,a) 

= d^^ ■ ^^'^ 

For a system of free particles we have ff^{t) = so that cr^*(i, a) — 0. If there are inter- 
particlc interactions, they will determine the effective stress tensor. 

Let us consider an arbitrary point fi{t). The multipole moments of the mass density p 
and momentum density pU and of the stress-like density pU^U"^ with respect to the point 
f]{t) are defined by setting (A^ > 0) 



N 

= E^i[ri?{T) - ri^^{t)]...[ril"{t) - ri^-{t)], 

i=l 

N 

n — m[?7(i)] — m — y^m^, 

i=l 

p^^--^^[fi{t)] = / d'a[a^^ -ri^^{t)]...[a^- -ri^-{tMt,a)U^{t,a) = 

V 
i=l 

N N 

n = Q p''[f]{t)] = J2 mmit) = I] < = ~ 0, 

1=1 i=l 

^ri...r„r.j^^^)| _ J d^a[a'' - 7]'' (t)]. ..[a'" - 7]'" {t)]p{t, a)^ {t, a)U' {t, a) = 



1=1 



The mass monopole is the conserved mass, while the momentum monopole is the total 
3-momentum, which vanishes in the rest frame. 

The point fjlt) is the center of mass if the mass dipole vanishes 



N 



i=l 



= 12^^[vl{t) - V''{t)] = ^ v{t) = q^r. (A8) 



The time derivative of the mass dipole is 
dm''[f]{t)] 



dt 



p'^iffit)] - mr]'{t) = fi+- mr^'it). (A9) 



When 77(f) = qnr-i from the vanishing of this time derivative we get the momentum-velocity 
relation for the center of mass 

P^[(inr\ = n\ = mq[^ [~ in the rest frame]. (AlO) 
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The mass quadrupole is 

N 



m'''[fj{t)] -^^miTjlitXit) - mrf{t)r]'{t) - {p'^{t)m'[f]{t)] + r)'{t)m'^[f}{t)]) , (All) 

i=l 

SO that the barycentric mass quadrupole and tensor of inertia are respectively 



N 



i=l 



u 



nri 



a.b 



N-1 



m 



'[qnr] = S''' J2 ^abPa ' Pb - 1''^% 



(A12) 



a, 6=1 



The antisymmetric part of the barycentric momentum dipole gives rise to the spin vector 
in the following way 

AT N 
i=l 1=1 



1 N-1 

S''=-e^'Y'[qnr]=Y.iPa^^,aT. 
^ a=l 



(A13) 



The multipolar expansions of the mass and momentum densities around the point fKt) 



are 



oo ^ri....rn\f^ an 



n=0 

oo 



p{t,a)U^{t,a) = Y: 



n=0 



pri....r„r^f^ Qn 



S'ia-rfit)). 



(A14) 



For the barycentric multipolar expansions we get 



Pit, a) = mS'ia - q^r) - ^inqur] - ^S^' E ^""[^-]) a^'^'^^ ~ ^""^^ + 



n=3 



-S^{a-qnr), 



Pit, a)U'(t, a) = - gk,) + [ie-"5« - q^r) + 



oo 



+ E 



n=2 



n! da^^...da^" 



(A15) 
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APPENDIX B: SYMMETRIC TRACE-FREE TENSORS. 



In the applications to gravitational radiation P5| , ^^ , ^ one does not use Cartesian ten- 
sorshut irreducible symmetric trace-free Cartesian tensors (STF tensors). While a Cartesian 
multipole tensor of rank I (like the rest-frame Dixon multipoles) on has 3' components, 
of which in general |(/ + + 2) are independent, a spherical multipole moment of order I 
has only 21 + 1 independent components. Even if spherical multipole moments are preferred 
in calculations of molecular interactions, spherical harmonics have various disadvantages in 
numerical calculations: for analytical and numerical calculations the Cartesian moments are 
often more convenient (see for instance Ref. for the case of the electrostatic potential). 
Therefore one prefers to use the irreducible Cartesian STF tensors ||49| , |50[| (with 21 + 1 inde- 
pendent components if of rank /), which are obtained by using Cartesian spherical (or solid) 
harmonic tensors in place of spherical harmonics. 

Given an Euclidean tensor Afc^...^^ on R^, one defines the completely symmetrized tensor 



Ski. 



ki 



A 



ir(l)---'«7r(/) ■ 



Then, the associated STF tensor is obtained by 



removing all traces ([//2] = largest integer < 1/2): 



[1/2] 



■^ki...ki ~ 2^ ^n0(kik2---0k2„-ik2„^k2n+i...ki)iiii...jnjny 
n=0 

n{2l-2n-l)\\ 



{I - 2n)\{2l - l)\\{2n)\\ 



(Bl) 



For instance (T^ 



abc) 



\STF 



(abc) 



ail) 
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